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Abstract 



We describe the absorption by the walls of a quantum electrodynamics cavity as a process during which the elementary 
O-i excitations (photons) of an internal mode of the cavity exit by tunneling through the cavity walls. We estimate by 
^— > classical methods the survival time of a photon inside the cavity and the quality factor of its mirrors. 
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1 Introduction 
> 

Atjsorption of a photon by the walls and/or escape of the 
imptons outside from a cavity is an intrinsically irrevers- 
il^ in time and non-unitary process. Its deep understand- 
ijlgj would require in principle an accurate description of 
properties of the environment (geometry and quality 
cG^^he mirrors in the cavity, possible escape towards the 
empty space outside the cavity and so on). However this 
4^^ription is too complicated to perform in the practice 
s^f^, commonly, one derives the effective master equation 
tl^ incorporates the losses inside the cavity thanks to an 
"educated guess". In this ad hoc description only one free 
parameter, the survival time of a single photon inside the 
cavity, characterizes the losses. The value of this parameter 
is fixed in last resort by experiments |1]. 

The aim of our work is to refine the usual description of a 
lossy cavity by 

(1) describing the absorption and/or escape by the bound- 
aries of a QED cavity as a tunneling process during which 
elementary excitations (photons) of an internal mode of 
the cavity exit the cavity by tunneling and 

(2) estimating by classical methods the survival time of a 



photon inside the cavity. 

The paper is organized as follows. 

Firstly (section [2| , we modelize the disappearance of the 
photons of the cavity as the tunnel process of a particle 
trapped between two square potential barriers (which rep- 
resent in our view the boundaries of the cavity). By doing 
so, we show that the classical picture according to which 
photons bounce back and forth inside the cavity is mislead- 
ing. Instead, the quantum state of light inside the cavity 
can be seen as a quasi-stationary state. This brings us 
to the next section (section [3]), where we refine this ele- 
mentary model by introducing a Quantum Jump Operator 
aimed at simulating the escape of the photons leaving the 
cavity. We show in the framework of this more sophistic- 
ated model that the survival time of a photon inside the 
cavity is in last resort defined by a single free parameter, 
which measures the coupling between a photon inside the 
cavity and a photon of same frequency leaving the cavity 
by tunneling through the barriers. This result is derived, 
in the weak coupling regime, along the lines of the so-called 
Friedrichs model which is very close in mind to the Wigner- 
Weisskopf treatment of exponential decay. 

An interesting by-product of our approach (section Q is 
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that it allows us to show in a straigthforward manner that 
coherent states couple to the environment without getting 
entangled to it. Usually this property is derived by solving 
the master equation that describes the losses in the cavity. 
In our case this property is seen as a direct expression of the 
(bosonic) indistinguishability of the photons in and outside 
the cavity. 

As, on the other side, coherent states are quasi-classical in 
the sense that they minimize Heisenberg-type uncertainty 
relations, it is consistent to neglect the quantum fluctu- 
ations in order to obtain a first estimation of the single 
free parameter that is not addressed by our quantum jump 
model, essentially the absorption rate of the cavity. In this 
coarse grained approach where the mode inside the cavity 
is treated as a classical oscillator, it is justified to resort to 
classical. Maxwell equations in order to estimate the losses. 
By doing so we find (section . [5] ) a good estimation of the 
lifetime of a photon inside the cavity, and obtain an indirect 
estimation of the quality factor of the mirrors used in the 
QED cavity |I7| developed by the group of Serge Haroche 
at the Ecole Normale Superieure de Paris [U |T6l [T71 [T^ . 



2 Tunneling of a particle trapped 
between two potential barriers 

A first possible model for cavity quantum electrodynam- 
ics is that of a quantum particle initially confined between 
potential barriers. Cavity losses are then modeled by the 
tunneling of the particle's wave function through the bar- 
riers. Initially we shall consider massive particles but later 
this approach will be generalized to particles of arbitrary 
mass, including zero mass particles, having in mind the 
application of our model to QED cavities. 

2.1 The Gamow model for nuclear decay 

The Gamow model |13)|14|fT5] describes nuclear disinteg- 
rations, namely, a-type disintegrations 
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Figure 1: Visualisation of the Gamow model. The central zone 
potential is purely strong force-induced, while it is purely electrostatic 
outside the strong force range. 



Both strong (nuclear) and electromagnetic (Coulomb) in- 
teractions were taken into account by Gamow, and the en- 
ergy Ea of the a particle was assumed in his model to 
be lower than the barrier height, so that it is in a "quasi- 
bound" state, but at the same time Gamow assumed that 
this energy was high enough so that the a particle can tun- 
nel through the Coulomb barrier towards the outside world. 
The situation is schematically represented in Fig. [T] The 
Gamow model is often [11][12| presented as follows: one 
takes for granted that the survival probability decreases 
exponentially, and is simply interested in computing the 
decay constant Vq. To this end, one writes it as a product 
of three factors 1121: 



Tg =p(a)/r 



(1) 



Helium nucleus 



where "p (a) is the probability that an a particle exists 
within the daughter nucleus but is independent from it 
(it is often taken to be equal to 1 since there is no easy 
way to evaluate it), / is the classical collision frequency 
between the a particle and the potential barrier created 
by the daughter nucleu^ and T is the tunneling probab- 
ility through the potential barrier. This can be summed 
up informally in the following way: the a particle oscil- 
lates back and forth inside the potential well, every now 
and then it hits the barrier and a small part of the wave 
function is transmitted through the barrier. This approach 
allows to derive the empirical Geiger-Nuttal law 

■^daughter 



log Vq = est 



(2) 



within a quantum mechanical framework. To this end, 
Gamow proposed a simple picture: within the mother nuc- 
leus, the a particle "sees" the potential created by the other 
part of the mother nucleus, namely the daughter nucleus. 



where is the energy of the a particle and ^daughter is 
the atomic number of the daughter nucleus. However, the 

^If one denotes by Ro the daughter nucleus's radius and Va the 
particle's speed, then we can write / — Va/ (2Ro). 



2 



Central region 




t n n 

Vo 

\ I I I _ 

—I — w —I I I + w 

Figure 2: The double potential barrier system. 



classical picture of a particle bouncing back and forth is 
misleading as we shall see and ought to be replaced by a 
purely quantum description. 

2.2 Results for a double potential barrier 

Inspired by the Gamow model, we study the so-called 
double potential barrier, sketched in Fig. [2] It features 
two potential barriers of height Vq situated between —l — w 
and — /, and I and / + w respectively. We are interested in 
the decay of a wave packet initially confined between the 
barriers. 

To solve the Schrodinger equation numerically, we followed 
[ lO] and defined the imaginary part and the real part of 
the wave function to exist at half-step time intervals from 
each other, in order to use the central difference method 
to compute the discretized derivatives. Remarks about the 
stability of the numerical algorithm can be found in |10| . 
Let us note that we had to add absorbing layers to avoid 
reflections. Following the perfectly matched layer (PML) 
method |32| 133). the Laplacian operator was thus analytic- 
ally continued in the complex plane [8]. 

In a first time we considered that the particle was pre- 
pared initially in a Gaussian wave packet (of energy de- 
noted i?free)- We let run different simulations for different 
values of the initial energy keeping constant Vq, the barrier 
height. Through numerical simulations we monitored the 
time evolution of the nonescape probability 

^'nonesc(i) = j dx {x , t)\^ (3) 

for different input energies |8]. Even though short-time 
behaviors are very different, namely, higher energy pack- 
ets decay much faster at short times, it turns out that at 
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Figure 3: Time evolution of the wave function square modulus for an 
initially Gaussian state. Gamow's bouncing picture is valid at short 
times, but that the wave function rapidly becomes a "quasistationary" 
wave that leaks through the barriers. 

long times, the decay is identical (and exponential) for all 
initial states. This means that the decay constant is in- 
dependent of the initial energy, and solely depends on the 
potential profile which the particle evolves in. At first sight 
this seems to be a hard blow dealt at the Gamow model, 
since we expect the collision frequency and the tunneling 
probability to depend on the "input" energy. On the con- 
trary, it appears that the Gamow model remains accurate 
if we focus on the specifications of the state at long times. 
The key step is to forget about the specifications of the 
initial state, and especially about the input energy, and to 
focus on the final, exponentially decaying state (which is 
the state the characteristics of which are accessible to the 
experimentalist who verifies the Geiger-Nuttal law). Ex- 
plicit numerical examples are given in [5]. 

Our numerical results are comforted by the theory of 
Gamow functions and complex poles, in particular by the 
leaky eigenmode approach |20j approach which consists in 
finding the system's normalized generalized eigenfunctions, 
which have poles in the complex k-plane. These poles cor- 
respond to the system's so-called Gamow functions, i.e. 
solutions to the time-independent Schrodinger equation 
for complex energy values. These complex energy eigen- 
values, the imaginary part of which is simply the corres- 
ponding Gamow functions' lifetime, are, up to dimensional 
factors, the square of the /c-values which are the normal- 
ized generalized eigenfunctions' poles. Using the tetracho- 
tomy algorithm (see [3 (sect. 3.4) or |8]) to find poles in 
the complex fc-plane, we were able to find several Gamow 
states and thus to observe departure from the seemingly 
"universal" behavior initially witnessed. In Fig. |4] we plot 
the time evolution of the two longest-lived Gamow states. 
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Figure 4: Time evolution of the wave function square modulus 
for the longest-lived (left) and next-to-longest-lived (right) Gamow 
states. One sees that the space-dependence of the wave function in 
the central region (that is, everywhere except in the outside region) 
is perfectly stationary and that the wave leaks through the barriers. 



We find that, accordingly to what is inferred in |20) . their 
space-dependence inside and between the barriers is time- 
invariant and that they simply "melt" with time. These 
Gamow states can thus be seen as leaky eigenmodes of the 
system for which the time-dependence can be factored out 
of a time-independent spatial envelope (see (Q). Let us 
note that they resemble the modes of an infinite potential 
well (the so-called harmonics), a fact that we shall exploit 
in the following section. 



Since Gamow's picture of a particle that bounces back and 
forth between the barriers proved rather irrelevant -the "fi- 
nal state" tends to be an (almost) standing wave which 
leaks through the barriers-, one might wonder why we find 
such a good agreement between the semi-classical model 
and the leaky modes approach. For instance, the estim- 
ation of the uranium decay constant obtained from ([T]) is 
Tg = 3.130 X 10~^^ s~^ j8], to compare with the value 
2.6722 Tpoie = 2.6722 x 10"^^ s"^ obtained from the com- 
plex poles treatment [20] and the experimentally observed, 
empirical, value Fexp = 2.5872 x 10~^^ s~^. Other examples 
of this kind are given in ref. |8|. In order to explain why 
it is so, we propose the following argument: at sufficiently 
long times, and in the central region (x G [— / — w,l + w]), 
we can write the wave function as 



■ip (x, t) = (f{x)e 2' 



(4) 



where, in first approximation, the longest lived state is, up 
to a global normalisation factor A^, proportional to the fun- 
damental mode of the potential well: (p{x) ~ Ncos{kx). 
It is then easily checked that 



dt J-l-TL 



dx IV' {x, t)\' 



jltldxiV'(x,t)r 



Making use of the probability conservation equation yields 
j{l + w,t) -j{-l-w,t) 



r 



where the probability current is given by 



j {x, t) 



h 

2mi 



^* ~dx ~ ^ 



(5) 

The wave function is purely outgoing at the interface 
between the potential barrier and the outside world, that 
is, at X = ±{1 + w), which allows us to write 



(x, t) = A- e-''"' e"^* at x 



-I — w. 



il) (x, t) = Aj^ e 



ikx -Ft 



at X = I + W 



and finally yields 



hk 



\A+\' + \A. 



2m{l + w) (|(^|2)„i,hin 



(6a) 
(6b) 



(7) 



where (•) within is the average taken over the central re- 
gion (see Fig. [2]) and k is equal to 27r/Ao, where Aq is the 
wavelength of the quasi-stationary mode inside the cav- 
ity, due to energy conservation. The first factor in this 
last equation is of particular interest, since it strongly 
resembles the bouncing frequency of the "semiclassical" 
Gamow model (see ([T])). As for the second factor, it is re- 
miniscent of Gamow's tunneling coefficient, since the latter 
links the intensity of the wave function at the outer end of 
the barrier to the intensity of the wave function at its inner 
end. 



3 Tunneling processes in terms of a 
quantum jump operator 

3.1 A quantum jump operator 

As we noted before, the Gamow states resemble the in- 
finite potential well eigenstates, which would motivate a 
perturbative treatment to the problem. However, the fact 
that the unperturbed system would feature infinite poten- 
tial regions is a major obstacle since one cannot easily "per- 
turbate oneself out of infinity". Accordingly, we turn to a 
different treatment where we couple the inside and out- 
side regions through an interaction Hamiltonian. This is 
a first step towards a second-quantization approach of the 
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problem which is a much more comfortable framework for 
photons. 

We make the approximation that the only relevant mode 
between the barriers is the one which resembles the longest- 
lived Gamow state. Between the barriers < /) it is 



equal to Vin {x, t) 



cos{kox) where ko = vr/ (2/) 



and hujQ is the (real part of the) energy of the longest-lived 
Gamow state. It is taken to be equal to zero inside the 
potential barriers and in the outside region (that is, for all 
l^l > /). For "escaping" wave functions outside the barri- 
ers we make the natural choice of the basis which consists 
of plane waves 'ip'^^(x,t) = e'C^^'-'^*)^ with A; > when 
X > w + I, and k < when x < —w — I. These func- 
tions are taken to be zero between and inside the barri- 
ers {—w — I < X < w + I). This choice is motivated by 
our interest in the decay of initially confined states. The 
Hamiltonian operator for the system finally reads 



(8) 



where i/Q°^°"* is a standard free-space quantum Hamilto- 
nian mutiplied by the indicator function for the adequate 
region. The quantum jump operator is defined through 



yjump 



dojik) [xiu{k))\rou:){i^- 

+X* {u;{k))\^,^){^P: 



out 



(9) 



At this level we do not make any specifications on the coup- 
ling function A, but it will be useful in the following to ask 
that it be sufficiently regular in the vicinity of ujq, and also 
that |A(a;o)|^/^^ is much smaller than ujq. As we shall 
see, the latter condition is typical of the Fermi golden rule 
regime, which is valid when the decay constant F is much 
smaller than the internal frequency lvq of the state. Our 
model belongs to the class of the so-called Friedrichs mod- 
els |22) in which a finite number {N) of discrete states is 
coupled to a continuum of decay products (here = 1). 

3.2 Direct sum and tensor product 

The natural structure of the system's Hilbert space Ti is 
given hy % = £2 (IS) = 'Hm ® l-Lout-, but here we turn to 
the different structure of the tensor product (this switch 
procedure from a direct sum structure to a tensor product 
structure is detailed in |22|): 



out • 



(10) 



We can further validate the procedure given in |22j with 
the following argument: for the direct sum structure one 
writes a state as 



{X I V'^'^) = C (^) llin (X) + (X) llout (X) 



di) 



(11) 



where Uin and llout are indicator functions for the in and 
out regions respectively. Thus the scalar product is given 

by 



out / • 

(12) 

Switching to the tensor product structure, we write 

(13) 

are normalized states. 



|V;«)=|^»)^|0out)+|0in)$3|VSt> 



where the vacuum states | Ojn/out) 
It is then readily found that one gets the same expression 



as (12) for the scalar product. We thus adopt this tensor 



product structure for the Hilbert space and write the state 
as follows: 



|V(t)) =a(t) |lin)® lOout) 

+ / dcoik) Pik,t) |0in)C5 \ujik)^J, 



(14) 



where a and /3 are complex amplitudes. Writing in this 
tensor product structure the Hamiltonian for the system 
which is equivalent to the one defined by ^ and ^ we 
get 

H =hwQ I lin)(lm| ®^out 

+ nin(» j Mk) h\u{k)^J{u{k)^^,\ 

Au:{k) [A {U {k)) I Oin) (lin I ® I W (fc)„^J (Oout I 

+A* (W {k)) I lin) (Oin I ® I Oout) {k),,, I] . 

(15) 

3.3 Time-dependent perturbation theory in 
the weak coupHng regime 

We shall now solve Schrodinger's equation for the Hamilto- 



nian defined by (15). Our treatment is very close in mind 



to the Wigner-Weisskopf perturbative approach. Since the 
state is initially confined inside the cavity, we choose 



VA; 



a(t = 0) 
/3 {k, t = 0) 



1, 
0. 



(16) 
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These initial conditions (16) yield the following equations: e.g., [23] for a more rigorous treatment) 



dt 



it)) 



dt 



fit) 



j dw (A;) 



\{^{k)) 



2 





gi(a;o-a;(fe))(t-t') j ^^l^ 



(17a) 



A(tJo) 



duj{k) 



\{uj{k)) 



h 

X{ujo) 



2 rt 





2 rt 



dt' / dL^(A;) ei(-o— (t') 



-TT 



A (g^o) 
h 



[ dt' 2tt 6 it - t') f it') 
Jo 



fit) 



(19) 



^17b) 



Let us note that up to this point we assume a bijective 
map from A; to w. This means that we only consider out- 
going modes to the right. We shall restablish the contribu- 
tion from the left-outgoing modes in the end. Now, since 
the "free" (non perturbed) Hamiltonian of the inside region 
reads Hi^ = huo \ lin)(lin |, we can reasonably assume, in 
the weak coupling regime, for which A is a small parameter 
and the decay rate is slow, that a is the product of an os- 
cillating exponential of frequency uq with a function / the 
variations of which are negligible on a time-scale of uJq : 



where the second step is approximately valid on energy 
conservation grounds. Under certain assumptions, f ap- 
pears as the solution of the first-order differential equa- 



tion (19). Together with (18) this motivates the following 



Wigner-Weisskopf exponential ansatz (with A = F -1- 2iwLS 
a complex- valued number): 



a(t) =e-^'^«*e-^^*. 



(20) 



a{t) 



fit). 



(18) 



This leads to (following a computation very similar to the 
one found in ref. p], sect. 2.5): 

1 r /|A (c^ {k))f\ 1 - e-i{'^W-'o)ie|At 
-A = 1 duj (k) - — J-^ -. . 



The first consequence of this assumption is that energy is 
conserved in first approximation, as can be seen by integ- 



rating (17a) considering that a{t) 
yields 



-iLUgt 



Collecting contributions from left- and right-outgoing 
modes, and in the A — ?■ O"*" limit the Sochocki-Plemelj the- 
orem j3] (sect. 8.1) allows to write, defining the transition 
f{t = 0) which rate 



4tt 



|A(u;o)r 



1/3 {k, t) |2 « I A/ (t = 0) (oj (k)) sin {{u - ojo) t/2) /h (uj - ojo)\ 



2 as well as the Lamb shift 



which can be considered to go to when t is large enough 
whenever (w— wq) > 27r/t which means that only the modes 
that vibrate at "nearly" the same frequency as the internal 



vp / do; (k) 



\^i^ik))r 

uj (k) — ujo 



(21) 



(22) 



mode get coupled to it by the Hamiltonian defined at ( 15 ). 



One can then reasonably assume that X(u> {k)) is constant 
and takes the values A (a;o) because anyhow the evolution 
will not couple the internal mode to external modes which 
would differ too much in energy. Coming back to this prop- 
erty at the end of the treatment it is easy to check that it 
is expressed by the weak coupling condition ojq ^T. This 



where vp denotes the Cauchy principal value of the sub- 
sequent integral, the following expressions ^ (sect. 3.5) for 
a (compare with (|19[)) and (3 [k, •): 



a(t) 
nk,t) 



-i(wo+WLs), 



r* 



(23a) 



o(k)t ^i^ik)) 1 

h 



_ g-i('^o+tJLS-t^(fc))ig-| 
UJ (k) - {UJQ + WLS 



rt 



allows us to write, developing the constraint (17b) (see. 



(23b) 
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3.4 Wavefront propagation 

It is straightforward to treat "relativistic" massless particles 
with our model (in the following we focus on bosonic 
particles) which allows us to derive in a simple fashion 
an expression for the spacetime dependence of outgoing 
wave packets in the Wigner-Weisskopf regim^ For mass- 
less particles these wave packets obey a linear dispersion 
relation uj (k) = ck. Let us now derive the expression for 
the "photonic wave function" in the outside region. It ad- 
mits a natural interpretation in terms of the photodetection 
probability function established by Glauber, as is shown in 
ref. |3], section 6.3, in the case of spontaneous emission 
of a two-level atom in free, three dimensional, space. We 
get, after a straightforward computation, (and modulo the 
assumption that all modes inside and outside are charac- 
terized by a same polarisation): 

V'right {x, t) = Uright {x) I dk(3 {k, t) e 



ikx 



V'left (X, t) = Ilieft (x) 



1 



dkf3{k,t) e 



-ikx 



(24a) 
(24b) 



The calculation (for details see [8]) yields 

^right/ieft(a::, t) = -^Aight/Mt{x)Q{ct =F x)A (wo) 

Thus we see that the wave function is nonvanishing inside 
the light cone only. At fixed time, it increases exponen- 
tially with increasing x or —x, depending on whether we 
look at what happens on the right or on the left of the 
inside region. The exponential growth is strongly reminis- 
cent of what is obtained in the framework of Gamow func- 
tions (for massive particles) |20) . as well as in the case of 
a non-quantum treatment of an open electromagnetic sys- 
tem (see [34j for more details). Although in the Friedrichs 
formalism, as far as we know, it is not common to give an 
explicit spacetime dependence for the decay products, this 
has been done in quantum optics jH |l6] and also in the 
leaky mode treatment of Maxwell's equations |31) . 

3.5 A Bohmian description of the trajectories 

Let us reconsider the right and left outgoing wave pack- 



where Uright and llieft are indicator functions for the right 
and left regions respectively. From the expression of f3 (see 



ets (27) and interpret them as the outgoing photonic 



(23)) we get (we extend the integration domains because 
the main contribution is centered around the pole) 



-(Aright (x, t) = llright (x) 



1 



^2tt J- 



+00 



dk 



Xick) 



ck-ujo + |A {ujo)[ 



^ik(x—ct) 



Ipleit (x, t) = llicft (x 



1 



27r 

ik{x+ct) 



-i(ujQt—kx) 



dk 



\{ck) 



(25a) 



ck- UJO + i^ |A(a;o)|' 



-i(Cjot+kx)-^\\[uio)?t 



(25b) 



wave functions. Interestingly, a simple Bohmian picture 
holds for these outgoing waves. Indeed the outgoing right 
(resp. left) wave function defined by ( [24| obey^5\I'/9t = 
±cd^ /dx, so that we can associate to the "wave function" 
^' a conserved density p{x,t) = \tjj{x,t)\'^ and a density 
current j {x,t) = zizcp{x,t). The Bohmian velocity reads 
then [19] t'Bohm(x,t) = j {x,t) /p{x,t) = ±c, and we re- 
cover the old Einsteinian picture of a photon moving at the 
speed of light in the outside region. Between the barriers, 
however, the particles are nearly at rest in the Bohmian 
picture, regardless of their mass, because the wave func- 
tion is nearly stationary in that region. This analysis con- 
firms our conclusion of sect. 2.2 in the bohmian picture. 



To compute these integrals, we use the Cauchy residue the- 
orem. We see that the pole Kq of the integrands is always 

1 



particles are at rest inside the barrers (cavity) and only 
begin to move once they leave it. 







UJQ 



^In the case of massive particles the result is essentially the same 
but we must then replace the velocity c by the group velocity of the 
particles at the energy fiLuo. Indeed, in the Fermi golden rule regime 
the dispersion relation is in good approximation linear for massive 
particles as well -with the de Broglie group velocity Vds = hko/m 
replacing c. This is so because for massive particles uj = hk^ / (2m) = 
h (fco + Skf' I {2m), and since in the Fermi golden rule regime tJo >> 
3uj « r, we have 5k << fco so that Slj = w — ljo — VdBSk. 



(26) 4 On the time evolution of coherent 
states 



In cavity QED experiments, photons are kept in a "box" 
(superconducting cavity) during "macroscopic" times (130 

■^This is also true in the case of massive particles, but then c should 
be replaced by the de Broglie velocity VdB = hko/m and j {x,t) is 
defined by ([5|. 
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milliseconds in the case of [IE]), after which they are ab- 
sorbed by the cavity walls. We shall now use the model 
defined in sect. [3]to study the behavior of coherent states in 
QED cavities, in particular their status of "classical pointer 
states". 



4.1 Decoherence and pointer states 

The decoherence program is aimed, among others, at ex- 
plaining why the macroscopic world does not (directly) 
obey quantum mechanical laws but instead behave classic- 
ally [25]. It is based on the insight that quantum systems 
are never isolated [27], but that their evolution is at least 
partly monitored by their environment with which they 
continuously interact (as is put in [1] (sect. 4), "the environ- 
ment is watching"). The decoherence program's goal is to 
build on this insight to explain how the "classical" proper- 
ties of the macroscopic world can be inferred from a purely 
quantum mechanical paradigm at the micropscopic level. 
In this approach, one can show |24| |25] how the interac- 
tion between the system and its environment destroys the 
coherence between a paticular set of quantum states, the 
so-called "pointer states", which correspond to the "pointer" 
positions of the macroscopic measurement apparatus, i.e., 
the different outcomes of a macroscopic measurement per- 
formed on the system [27]. The so-called quantum pointer 
states are "quasi-classical" states of a quantum system. 
They are singled out by their robustness to the interac- 
tion between the system and its environment. Since deco- 
herence is, as is well-known, the corollary of entanglement 
[2T| |28], pointer states of the system can be defined as the 
states which become minimally entangled with the envir- 
onment in the course of their evolution [26j. In our model, 
the inside region plays the role of the system while the 
outside region mimicks the environment. 

4.2 Decoherence-free behavior as a con- 
sequence of bosonic undistinguishability 

Let us write symbolically 



(28) 



We know from sect. 3.3 that an initially confined state 
I lin)® |Oout) evolves into a superposition of the type 

a (t) I lin)«) I Oout) +Ht)\ Oin)® I lout). 



Thus, we can write the following time evolution for n (un- 
distinguishable) bosons which we label with the index k: 

n 

® (|lin)fc«) |Oout)fc) 



fc=l 



(a (t) I lm)fe^ I Oout)fc + /3 (t) I Oin)fe«) I lout>fc) 



k=l 



(t) P^-^ (t) 



j=0 



j\ {n-j)l 



\3ir 



out / ' 



where the normalized states 



[n-J, 



out/ 



(29) 



are 



defined by symmetrization over all possible states having 
j particles in and {n — j) particles out. The corresponding 
nonnormalized states read 

[|jin)«>| (ra-j)out)]non 

= ((I lin>l«> I Oout>l)® • • • «'(| lin>j® I Oout)i)) 

(g) ((|0in)j+lO| lout)i+l)® • • • O(|0in)„(» I lout)n)) 

-|- all other combinations with j in and (n — j) out. 

There are (nl) / {j\ (n — such combinations, which 
means that 



|jin)®| (n- j)out) 



j\ {n-jy. 



It is straightforward to see that 



out / J non ' 

(30) 



(a {t) I lin)fc«' I Oout>fc + /3 (0 I Om)fc^ I lout)* 



fc=l 



^a^-(t)/3"-^-(t) [|iin)®|(n-i)„„,)U 
i=o 



which along with (30) proves (29), which could be labeled 



a "quantum binomial" law. In particular, if at time t = 
the cavity is prepared in a coherent state and the exterior 
is in its vacuum state, then the evolution generates a tensor 
product of an interior and an exterior coherent states: 



+ 00 



|^(t = 0))=e'^ 



n=0 



l^'in)'^ I Oout) 



+00 n 



i^(o)=e-^H*)i^j:5: 



> e 2 



\m\' 



n=0 m=0 

\n—m 



y^(n — m)! 



(n — m)^ 



(31) 
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where we used the identity |a + |/3 
at time t can be rewritten as 



1. The state 



|^(t))=e-^l°WI^V ^^"ff)' |U 
,^ V A;! 



k=0 

+00 



I e 2 



M!,„„p|,(iml:! ,,„„„.,. (32, 



m=0 



This estabhshes that coherent states of the cavity interact 
with the exterior without getting entangled with it. They 
can thus be considered as "classical pointers" according to 
the criterion for classicality derived by Zurek in [26J. A 
common argument for the classicality of coherent states is 
based on the Lindblad master equation, which is a widely 
used dissipative (i.e. nonunitary) equation which mod- 
els the coherence loss which a quantum system undergoes 
when it interacts with its environment. It governs the stud- 
ied system's reduced density matrix, and, in the case of a 
QED cavity (at zero temperature), it reads [Ij (sect. 7.5.1) 



dt 



pin (t) 



+ 



1 

ih L 

r 



Hin, Pin (t) 

2apui (t) a) — a^apin (t) — pin (t) a^a . (33) 



In [21] we proposed an alternate derivation of the Lindblad 
master equation. Our proof is based on the fact that any 
damped coherent state of the form 



I6n(t)) 



i"(t)r 



+00 

E 

m=0 



rriw 



(34) 



with a given by (23) and H = huja^a, obeys the Lindblad 



equation (33). Since coherent states constitute a basis of 
the Hilbert space (rigorously, they form an overcomplete 
basis), we find that, by virtue of the linearity of the master 
equation, any state is a solution of the Lindblad master 
equation. 

4.3 Damped coherent states and the classical 
limit 

It is sometimes argued that the Ehrenfest theorem estab- 
lishes that quantum average values behave like classical 
quantities. This is wrong because, even if one can de- 
rive in a variety of situations an equation that looks clas- 
sical at first sight, of the type md^(r) /dt^ = (F), this 
equation is not equivalent to its classical counterpart in 



general. For instance in the case where the force derives 
from a potential, F = -VV (r, t) and (F) = - {VV (r, t)), 
but, generically, (VV (r, t)) 7^ ((r) , t) and it is only in 
situations where we can neglect quantum fluctuations that 
classical mechanics can be inferred from quantum mech- 
anics, which is in a sense a tautology. Damped coherent 
states however constitute an exception in the sense that 
whichever the state of a damped oscillator could be, we 
have m (P (x) / dt"^ = F {{x) (t) ,d{x)/dt). This is easily es- 



tablished by estimating thanks to the master equation (33) 



the temporal evolution of (a). Elementary computations 
establish that 



4^= -w--|{<»(*). 



(35) 



After integration over time we get 

(a) (t) = (a) (t = 0)e(-''^-?)*. 

Making use of the fact that (x) (t) = (a) (t)-|-(d)* (t), we get 
finally that (x) (t) = ^dcos {yit -|- 0o) exp ((— r/2) t), which 
is the most general solution of the dynamical equation 
of a classical-one dimensional-damped oscillator (that is, 
d^x/dt^ + rdx/dt-FcIi^x = 0, with ufl = From 
this point of view, it is justified to describe the average 
behavior of a damped oscillator trapped in a QED cavity 
as a classical oscillator, and thus to resort to the Maxwell 
equations in order to study its properties, as we shall do 
in the next section. It is by the way particularly relevant 
to approximate the behavior of the quantum oscillator by 
its average behavior in the case of coherent states because 
they minimize Heisenberg uncertainties, which emphasizes 
their status of "quasi-classical states". 



5 Quasimodal analysis of a QED cavity 

In this section we study the open QED cavity described in 
|17) by searching for its eigenmodes and complex eigenfre- 
quency, using a Finite Element Method (FEM). This clas- 
sical electrodynamic approach allow us to derive a num- 
ber of features observed experimentally by Haroche and 
coworkers (|18|). 

5.1 Searching resonances of open resonators 
using perfectly matched layers 

The major difficulty in the treatment of open problems in a 
numerical scheme based on a finite computational window 
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is to deal with infinity issues. Since their introduction 
by Berenger in [35] for the time dependent Maxwell's 
equations, Perfectly Matched Layers (PMLs) have become 
a widely used technique in computational physics. The 
idea is to enclose the area of interest by surrounded layers 
which are absorbing and perfectly reflectionless. These 
absorbing boundary conditions can be understood in the 
global framework of transformation optics (L38J)- The 
principle of the technique is to perform a geometrical 
transformation (here a complex stretch of coordinates), 
leading to equivalent material properties ([5l [36l 137]). 

The spectral problem we are dealing with consists in find- 
ing the solutions of source free Maxwell's equations, i.e. 
finding complex eigenvalues = {con/c)'^ and non zero 
eigenvectors such that : 

M{En) := V X (^-1 • V X = • En- (36) 

where e and /i are the relative dielectric permittivity and 
magnetic permeability tensors describing the electromag- 
netic properties of the system (cavity+external world). 
For Hermitian open problems, the generalized spectrum of 
Maxwell's operator M is real and composed of two parts 
: the discrete spectrum with trapped modes exponentially 
decreasing at infinity, and the continuous spectrum with 
radiation modes oscillating at infinity. In addition, another 
type of solution is present and very useful to characterize 
the spectral properties of unbounded structures : the so- 
called leaky modes (also termed quasimodes, quasi normal 
modes ([5l[39]) or quasi guided modes (|40j) in the literat- 
ure). These eigenmodes with complex associated frequency 
are an intrinsic feature of open waveguides. 
PMLs have proven to be a very convenient tool to com- 
pute leaky modes in different configurations (|4H I42j). In- 
deed theyefficiently mimic the infinite space provided a 
suitable choice of their parameters. The introduction of 
infinite PMLs rotates the continuous spectrum in the com- 
plex plane (since the operator involved in the problem is 
now a non self-adjoint extension of the original self-adjoint 
operator). The effect is not only to turn the continuous 
spectrum into complex values but it also unveils the leaky 
modes in the region swept by the rotation of this essential 
spectrum ( |9] ) . It is important to note that leaky modes do 
not depend on the choice of a particular complex stretch- 
ing : adding the PMLs is only a way to discover them. 
Finally, in order to apply the FEM, the PMLs have to be 
truncated at finite distance which results in an operator 



having only point spectrum with approximate radiation 
modes (also termed as PML modes or Berenger modes) 
due to the discretization of the continuous spectrum by fi- 
nite PMLs (|43)). Actually, the PML technique is analog 
to the theory of analytic dilation ( |44| I45| ) developed in 
the 1970s for Schrodinger Hamiltonians. 
In the sequel, eigenvalues are denoted w„ = cj^ -|- iuj'^. The 
real part is the resonant angular frequency u'^ = 271 fn and 
the imaginary part is the damping coefficient, which is 
related to the lifetime r„ of the photon in the cavity by 
cj^ = 27r/r„. The quality factor associated to a resonance 
is defined by Qn = oj'n/ (2^^)- 



5.2 Numerical results 
5.2.1 Numerical set up 

A detailed description of the cavity can be found in |17) . 
It is composed of two mirrors of diameter D = 50 mm 
facing each other. The distance between their apexes is 
L = 27.57 mm, and their surface is toroidal with radii of 
curvature r = 39.4 mm in the Oxz plane and R = 40.6 mm 
in the Oyz plane. The mirrors are coated with a thick 
layer of superconducting niobium. We take advantage of 
the symmetries of the problem and thus model only one 
eighth of the cavity to save memory and computation 
time. By setting a well chosen combination of Neumann 
and Dirichlet boundary conditions on the cutting planes, 
we can select the modes with desired symmetries. The 



eigenproblem defined by (36) is then solved by the FEM, 



using cartesian PMLs terminated by homogeneous Neu- 
mann boundary conditions to truncate the infinite space. 
The computational cell is meshed using second order 
edge elements, with a maximum size of an element set 
to y / (|lHee|A^), where = 3.68 cm is the approximate 
resonant wavelength of the cavity, and is an integer 
(A^ = 7 for the domain inside of the cavity, N = 7 for 
the domain outside the cavity, N = 5 for the PMLs and 
A'^ = A^m for the mirror surfaces). The final algebraic 
system is solved using a direct solver (PARDISO). 

In a first approximation, both mirrors are assumed to be 
lossless, thus we set Perfect Electric Conductor boundary 
conditions (PEC, i.e. Dirichlet boundary conditions). In 
a second step, in order to account for losses, absorption is 
considered through a Surface Impedance Boundary Condi- 
tion (SIBC) on the boundaries of both mirrors. 
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5.2.2 A first approach: Perfect mirrors 



Let us consider first a lossless model, i.e. perfect mirrors. 
Among many others, we find only two eigenfrequencies 
with exceptionally low imaginary parts (or high associ- 
ated lifetimes). Their real parts correspond to frequen- 
cies fi and /2 located around 51.085 GHz (to be com- 
pared to 51.099 GHz found experimentally in |17|). with 
r = /i — /2 ~ 1.2 MHz(same as in [17J). According to the 
well known Fabry-Perrot interferometry principle, a modi- 
fication of the distance between the mirrors would allow to 
retrieve the experimental frequency, as shown in the fol- 
lowing section. The small discrepancy between the two 
complex eigenfrequencies is due to the removed degener- 
acy induced by the slight cylindrical symmetry breaking 
of cavity (r ^ R). The mode labelled 1 with higher fre- 
quency is x-polarized while mode 2 with lowest frequency 
is 7/-polarized, and they both have nine antinodes along Oz 
(modes TEMgoo, see the field maps on Fig. [s]). 

We study the convergence of the eigenfrequencies as a func- 
tion of the mesh refinement. A satisfying convergence is 



obtained with increasing (Fig. 6a) for the real part 



of the frequencies, while all that can be said about life- 



times (Fig. 6b) is that they have an order of magnitude 
of a few seconds. This non-convergence is due to the fact 
that the resonances features are in that case dominated by 
the geometry of the mirrors surface, which is approxim- 
ated by tetrahedral elements, thus not exactly reproducing 
the smooth curvature of the toroidal surfaces. To be more 
specific, Haroche et al. fabricated two cavities Mi and 
M2. They measured a fabricated maximum peak to val- 
ley deviation from the ideal toroidal mirror shape of only 
300 nm. We can easily estimate a numerical maximum 
peak to valley deviation of our virtual tessellated mirrors. 
The 3D mesh is made of tetrahedra, so the surfaces of in- 
terest of each mirror are discretized into a set of triangles 
of side lengths set to a maximal value of \^ / {\yKte\N). 
Let us consider an equilateral triangle T of side length 
y / (|9^ee|A^) whose vertices belong to a sphere of radius r. 
When approximating the corresponding piece of sphere by 
this triangle T, the maximum deviation from the original 
shape is the minimum distance between the barycenter of 
T and the sphere. For Nm = 22, this characteristic length 
is approximately 300 nm, which corresponds to the exper- 
imental value. Our numerical mirrors are globally closer 
from the exact toroidal shape than the experimental ones. 
However, numerically taking into account the remaining 
local surface roughness of typically 10 nm is not only nu- 
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Figure 6: Convergence of the eigenfrequencies with mesh parameter 
for modes 1 and 2: (a) resonant frequency /, (b) lifetime r 



merically unthinkable but also irrelevant (details of Ar over 
half a million). Finally, the lossless PEC model represents 
a cavity so resonant that even a tiny change in a fine (set- 
ting Nm to 60 leads to a maximum deviation from the 
exact shape of only 40nm, i.e. Ar/150000!) discretiza- 
tion leads to a significantly different lifetime. In spite of 
the non-convergence of this first model, the order of mag- 
nitude of the lifetime is far greater than lifetimes found 
experimentally (seconds instead of 100 ms in [17]). A first 
conclusion can be drawn: Assuming a perfect geometry 
(perfect mirrors with ideal shape and alignment), losses 
are not dominated by radiation loss. We now try to give 
an estimate of the effective resistive losses. From a nu- 
merical point of view, including some friction mechanism 
is expected to benefit to the model convergence. 



5.2.3 A more realistic approach: Lossy mirrors 

To that extent, we now apply Surface Impedance Bound- 
ary Conditions (SIBC) = Xs+ lY^ on the mirrors faces. 
London penetration depth for niobium Ll (independent 
of frequency) is set to a typical value of 0.1 ^m [6]. In 
the framework of a two fluids model, at low temperature, 
the imaginary (inductive) part of the impedance can be 
approximated by Kj = w/io-f'L = 6.4 jsQ. As for the real 
(resistive) part of the impedance, it is extremely difficult 
to measure and greatly depends on numerous experimental 
conditions. Therefore, our only option is to estimate this 
parameter to obtain lifetimes of the same order of mag- 
nitude as measured in jT^. Finally, we adjusted the length 
of the cavity to find resonant frequencies closer to those 
measured in [17]. 

With the updated value of L = 27.562 mm (instead of 
27.57 mm in llTj), the system exhibits two resonant fre- 
quencies /i = 51.0984 GHz and /2 = 51.0997 GHz (r = 
/i - /2 1.29 MHz). As in the PEC model, the con- 
vergence as function of mesh refinement is reached (see 
Fig. [7^ for A^^ = 40. 



Moreover, the discrepancy between lifetimes when refin- 
ing the mesh size on the mirror is greatly reduced com- 
pared to the PEC case (see Fig. 7b). With Xg = 1 /ifi, we 
obtain lifetimes of ~ 100 ms for both modes, which corres- 
ponds to the average lifetime found in [TT] for both cavities: 
112 ± 4ms (LF mode) and 87 ± 10ms (HF mode) for Mi, 
74 ± 6 ms (LF mode) and 130 ± 4ms (HF mode) for M2. 
This value of should be seen as an upper bound of the 
resistive phenomenon. It includes all loss processes beyond 
radiation loss: roughness of the mirrors, superconductor 
imperfections. . . 

Finally, it is indeed much easier to "align the cavity numer- 
ically" than experimentally! The fact that the discrepancy 
between lifetimes corresponding to the two cavity modes is 
greater experimentally than numerically tends to indicate 
a residual tilt (3 Euler angles) between the two mirrors. 



6 Conclusions and open questions. 

Many ideas exploited by us in the present paper are already 
present in classical quantum optics textbooks. For instance 
it is well known that, when an oscillator is linearly coupled 
to an arbitrary number of oscillators, coherent states factor 
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Figure 7: Convergence of the eigenfrequencies with mesh parameter 
TVin for modes 1 and 2: (a) resonant frequency /, (b) hfetime r 



out during the temporal evolution. In other words, coher- 
ent states are the quasi-classical pointer states of the prob- 
lem. However our approach is slightly different because, 
among others, we treat the problem in a first quantization 
approach. This has as a consequence that it could be gen- 
eralised to slow, massive particles, cold bosonic atoms for 
instance. It also sheds a new light onto certain "weird" 
features of quantum optics like the so-called Hong-Ou- 
Mandel two-photon interference effect |29j which is usually 
presented in terms of annihilation and creation operators. 
Roughly speaking, the effect is the following: whenever 
two indistinguishable photons are sent along two arms of 
a balanced (50-50) beamsplitter (schematized by the fol- 
lowing diagram: CutA)' they interfere constructively 
in the outgoing arms so that both photons belong to the 
same ongoing arm. Schematically, the beamsplitter acts as 



follows: it sends an ingoing state 



-iin im\ 



onto the state 



I 2^*, 0^"H I Oa *> 2b *)) /\/2', where the states | i^) (Ns)) 
i = 1,2) are generated by letting act the creation operat- 

t fj 



ors a 



''A {'^'b) ^ times for the A (B) mode respectively. 

In our approach, the initial state is the properly symmet- 
rized two boson state 



1 



1 ™ 
^A 



,0 



b) 



As described in |1J , the passage through a beamsplitter can 
be modeled by the coherent exhange of quanta between 
the two modes A and B. In the present case, its effect 
resorts, in the case of a balanced beamsplitter, for which 
the interaction time is the fourth of the Rabi period, to the 
transformation law |21] 

I n, OS) ^ (I ir , or )+ 1 or, ir » /V2, 
I o^A, m ^ (I ir , or)- 1 or, ir)) /V2. 

The passage through the beamsplitter is thus expressed, in 
this first quantization scheme, by the transformation law 
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The effect is obviously exactly the same as the Hong- 
Ou-Mandel effect, but no increase of entanglement ap- 
pears during the beamsplitting process: the inital and fi- 
nal states are both maximally entangled Bell states in this 
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picture, contrarily to what happens in the second quant- 
ization formulation where the initial state is a factorizable 
state. Although this program lies out of the scope of the 
present work, it would be highly interesting to figure out to 
which extent our alternative, first quantized, formulation 
of quantum optics could be pursued. 

Another non-standard feature of our treatment is that we 
never introduce the concept of quantum jump. 

This contrasts for instance with a popular approach for de- 



riving the master equation of a lossy QED cavity ( 33 ) , the 



so-called Monte-Carlo approach [T| [SO] according to which 
(here we consider the temperature of the "environment", 
that is to say of the modes outside the cavity, to be equal 
to zero), after an infinitesimally small time increase dt, 
either the mode inside the cavity loses one photon (with 
probability Fdt), or it remains unaffected by the absorp- 
tion (with probability 1 — Tdt). This random process can 
be simulated as a Monte-Carlo process and it leads to the 



derivation of the master equation (33) after averaging on 



the two possible outcomes |21j (either a detector 'clicks" 
and reveals the presence of an outgoing photon or it does 
not click). Certain features of this model are somewhat 
puzzling. For instance, in the case where the initial state 
of the quasi-stationary mode inside the cavity is a coherent 
state, it disentangles from the outside world throughout its 
temporal evolution, and it is difficult to figure out how the 
measurement of a photon outside the cavity would affect 
the state of the mode inside the cavity, even at the mere 
informational level, because the corresponding modes are 
totally decorrelated. 

In contrast, our approach makes it possible to provide a co- 
herent, unitary, description of the evolution of the full sys- 
tem (modes inside AND outside the cavity). In our view, 
decoherence (which manifests itself through the Lindblad 



equation (33)) is simply seen to be a corollary of entan- 
glement and during the whole process, the state of the 
full system (cavity+external world outside the boundar- 
ies) remains pure and evolves unitarily. Coming back to 
the famous Schrodinger cat problem, where a cat trapped 
in a box is prepared into a superposition of a living and of 
a dead cat, our approach makes it possible to give a com- 
plete description of the dead-living cat AND of its direct 
environment, inside the box. In this view, it is only when 
we open the box that the superposition breaks down. In 
a sense, our approach thus makes it possible to provide a 
complete description of what happens inside the box. It is 
only when we open the box that the superposition possibly 



breaks down. 

Besides these foundational issues, the classical techniques 
that we used for estimating the survival time of a photon 
inside the cavity (sect.jsj) confirm that classical electromag- 
netism is sufficient to obtain the main features of the losses 
inside the cavity. For instance our numerical simulations 
also predict that the two longest-lived modes inside the 
cavity are characterized by orthogonal polarizations and 
possess 9 nodes, and they provide a quantitatively correct 
estimation of their frequencies as well as of their frequency 
separation. 



As we have explained in sect. 5.2 , our estimations are more 
qualitative than quantitative for what concerns the life 
times of these modes but at least we derive a correct or- 
der of magnitude in that case. Even then, our numerical 
treatment delivers some information relatively to the losses 
of the cavity, information that would be difficult to gather 
by other methods. For instance, our analysis suggests that 
the main source of losses is not the escape of the photons 
to the free space which surrounds the cavity but is rather 
due to the absorption by the mirrors. Similarly, it reveals 
that the rougness is not a weakness of the cavity regarding 
losses but that the alignment of the mirrors is a crucial 
parameter. 



Finally, our results outlined in sections sects. 2 , 3.4| and [3751 
suggest to abandon the naive classical picture of the photon 
that is represented as a pointlike particle which oscillates 
back and forth between the mirrors and to replace it by 
the more appropriate picture where the photon inside the 
cavity is associated to a quasi-stationary standing wave 
trapped between the two mirrors. 
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